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The probability distribution of finite-time Lyapunov exponents provides an important charac- 
terization of dynamical attractors. We study such distributions for strange nonchaotic attractors 
(SNAs) created through several different mechanisms in quasiperiodically forced nonlinear dynam- 
ical systems. Statistical properties of the distributions such as the variance and the skewness also 
distinguish between SNAs formed by different bifurcation routes. 
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I. INTRODUCTION 



The characterization of attractors in nonlinear dynamical systems is a problem that has seen considerable progress 
in the past decade Beyond classification as simple or strange, the calculation of fractal dimensions and the 

description of the measure in term of a multifractal spectrum of singularities has become an important means of 
describing the structure of dynamical attractors [Q . For hyperbolic attractors there are rigorous results connecting 
the multifractal structure, through the thermodynamic formalism, with dynamical information as embodied in finite- 
time Lyapunov exponents . 

Grebogi et al. [Q first described dynamical systems wherein the attractors that result are fractal, but the dynamics 
is not chaotic, in that the largest Lyapunov exponent is not greater than zero. These strange nonchaotic attractors 
(SNAs) are generic in quasiperiodically forced systems. Subsequently, considerable effort has been directed toward 
the characterization and study of SNAs [|5|-|T]| , which have also been observed experimentally p2|-p[ . A potential use 
of SNAs is in the area of secure communications, and recent applications |r^ , [l6[| exploit the ease of synchronization 
of such systems. 

The strange nonchaotic state is only one of the possible dynamical states realized in quasiperiodically driven systems; 
periodic, quasiperiodic and chaotic attractors can also be obtained as parameters are varied. SNAs are typically found 
for parameter values very close to the boundaries of the chaotic regions, and the different bifurcation mechanisms 
through which they are created is a problem of interest. There are a number of routes or scenarios for the creation of 
SNAs, some of which can be correlated with bifurcations. These include 

(i) the Heagy-Hammel (HH) [|| mechanism involving a collision between a period-doubled torus and its unstable 
parent, 

(ii) the blowout bifurcation route and 

(iii) intermittency Jl7t , when as a function of driving parameter a chaotic strange attractor disappears and is 
eventually replaced by a torus through an analogue of the saddle-node bifurcation. 

The signature of these bifurcations in terms of the behaviour of the largest Lyapunov exponent has been discussed 
in detail The blowout bifurcation mechanism requires that the quasiperiodic torus of a system with an 

invariant subspace losses its transverse stability as a parameter changes across the transition and lead to the birth of 
an SNA. In this process the transverse Lyapunov exponent becomes positive while the nontrivial Lyapunov exponent 
for the whole system remains negative. In the HH mechanism Q, as system parameters are varied, the period-doubled 
torus gets progressively more wrinkled and collides with a parent unstable torus; this scenario is like an attractor 
merging crisis p9| . The distinctive signature of the intermittency route to SNA is a sharp change in the Lyapunov 
exponent which shows large variance and scaling behavior at the bifurcation. 

A general mechanism that is frequently observed but for which there is no well-identified bifurcation is the so-called 
fractalization route ||^, whereby a smooth torus gets increasingly wrinkled and transforms into a SNA without any 
interaction with a nearby unstable periodic orbit (in contrast to HH). This is probably the most common route to 
SNA in a number of maps and flows 

The present paper addresses the issue of distinguishing among SNAs formed by different routes through the use of 
finite-time (or local) Lyapunov exponents. We show that the morphologies of different SNAs differ in crucial ways, 
particularly for intermittent SNAs 1 17 1^ . This is seen most dramatically in the characteristic distributions of local 
Lyapunov exponents and the statistical properties of the distributions such as the variance and the skewness. 

In Sec. II, we briefly introduce the dynamical systems that are studied here. Results are discussed in Sec III. This 
is followed by a summary in Sec. IV. 



II. DYNAMICAL SYSTEMS 



Several quasiperiodically driven systems — both maps and flows p|-pl| — have been shown to have SNAs. We consider 
the quasiperiodically forced logistic map wherein three of the routes to SNAs can be observed. This system is 
defined by the equations 

Xn+i = a[l + ecos(27r0„)] x„ (1 - a;„), 

</>«+! =0n+t^ (modi), (1) 

where x £ li^ , (I) & , uj — (a/S — l)/2 is the irrational driving frequency, and e represents the forcing amplitude. 
The largest nontrivial Lyapunov exponent A is defined through 



2 



1 ^ 

AAr(x„) = — y^ln \dxn+j /dxn+]-i\, (2) 
A = lim AAr(x„), (3) 

Aat being the local or A^-step Lyapunov exponent. Note that Ajv depends on the initial condition a;„ while A, of 
course, does not. 

A region of the phase-diagram of the forced logistic map, Eq. (|^) in the a - t' plane ( where e' — e/(4/a — 1)) 
is shown in Figure 1. The rescaling of the forcing amplitude is done for convenience and the region of Figure 1 
corresponds to the period-3 window of the unforced logistic map. The symbols P, S and C correspond to periodic, 
strange nonchaotic, and chaotic behavior of the system. The dashed line is the locus of the period doubling bifurcation 
from period-3 to period-6. 

Along the left edge separating the periodic window from the region Ci, there appear to be no SNAs. On the right 
edge marked C'^, fractalized SNAs are obtained. The intermittent SNAs are found in the region marked I. Within 
this window, it was difficult to locate the HH mechanism for the formation of SNAs which is known to operate in this 
system for a number of different parameter values [p|jl7| . 

We also study the blowout bifurcation route to SNA, which occurs in the mapping [nij 



Xn+i — [a cos(27r(/)„) + h] sin(27ra;„) 

(t>n+i ^ (f)„+uj (mod 1), (4) 

(a and b are parameters and uj = {^/b — l)/2) which is bounded and has no other stable attractors other than the 
invariant subspace (x — 0, cf)). As parameters changes across the critical values, the dynamics of x leads to a strange 
nonchaotic attractor JTl] ! . Our interest is in contrasting this mechanism for SNA formation which is also accompanied 
by on-off intermittency, with the intermittent SNA p7| . 



III. RESULTS 



Although Aat depends on initial conditions, the probability density, defined through 

P{N, X)dX = Probability that lies between A and A + dX, (5) 

does not. This distribution can be obtained by taking an (infinitely) long, ergodic trajectory, and dividing it in 
segments of length N, from which the local Lyapunov exponent can be calculated through Eq. (2). 

For chaotic motion it has been argued [^,^ that since the local Lyapunov exponents can be treated as independent 
random fluctuations, the central limit theorem is valid, leading to a normal distribution for Aat, 

P(iV, A) « , ^ exp[-iVG"(A)(A - A)V2] (6) 
^2^iVG"(A) 

with the function G, the spectrum of effective Lyapunov exponents Q , being appropriately defined Q . 

These expectations are not always satisfied since there can be important correlations in the dynamics. We have 
recently described the characteristic distributions for finite time Lyapunov exponents in low dimensional chaotic 
systems where there are significant departures from central-limit behaviour [po[ . 

For SNAs there are additional complications. Although A is by definition nonpositive, P{N, A) can have a significant 
contribution from A > 0: for some of the time, these systems behave chaotically because of the fractal structure of 
the attractors 0. In the limit of large N, the contribution from positive A decreases and the density collapses to a 
(^-function, limAr^oo P{N, A) — > (5(A - A). 

Shown in Figs. 2 and 3 are local Lyapunov distributions for the four routes to SNAs (parameters are specified in 
the caption), for short (N = 50) and long but finite {N = 1000) times. The fractalized and HH SNAs both show a 
gradual approach to the normal distribution (a gaussian is fit to the data in Fig. 3), while the blowout SNAs, and 
more spectacularly, the intermittent SNAs show a distinctive departure from the gaussian distribution. 

The intermittent SNA is morphologically and dynamically very different from the other SNAs, and the shape of the 
characteristic distribution, P(N,X) is a combination of a gaussian and an exponential ||2^. In contrast to the other 
SNAs, the distribution is asymmetric, and the large A tail decays very slowly. 

It appears that the intermittent SNA is in a distinct universality class [ pO[ and a number of quantitative measures 
can be devised in order to show this distinction. Consider, for example, the fraction of exponents lying above A = 0, 
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F+{N) = / F(iV,A) dX, 
Jo 



(7) 



F+{N) vs N for the different SNAs are shown in Fig. 4. Except for the intermittent SNA, for which i^+(A^) ~ N~^, 
this quantity decays exponentiaUy, F^{N) ~ exp(— 7iV), with the exponents /3 and 7 depending strongly on the 
parameters of the system. We found that the values of the exponents are (3 = 0.72 for the intermittent SNA and 
7 = 0.02,0.007 and 0.042 respectively for HH, fractalized, and blowout SNAs. Similarly, other statistical properties 
of the these distributions can be studied. We calculate the first two moments about the arithmetic mean of all 
distributions and obtain the variance, 

/oo 
{X- Af P{N,X) dX, (8) 
-00 

and the coefficient of skewness, namely 

/oo 
{X-Af P{N,X)dX/{a^), (9) 
-00 

which are shown in Fig. 5(a) and Fig. 5(b) respectively. Generally for all types of SNAs, the variance of P{N,X) 
decreases as a power of N, ~ where the exponent S is different for each SNA. The variance for intermittent 

SNAs decreases very slowly compared to other SNAs, and our numerical results for the exponents, for the exam- 
ples shown here are = 0.97,1.71,1.63, and 1.7 for intermittent, fractalized, HH, and blowout SNAs. The degree of 
asymmetry in the distribution is quantified by the significantly larger skewness 71 (see Fig. 5(b)). 



IV. SUMMARY 



In the present paper we have studied the dynamical structure of strange nonchaotic attractors formed by different 
bifurcation mechanisms in quasiperiodically driven systems, by examining the distribution of finite-time Lyapunov 
exponents. Although the Lyapunov exponent is negative on a SNA, over short times, nearby trajectories can separate 
from one another since the attractor is strange: this corresponds to a local positive Lyapunov exponent. The manner 
in which this distribution changes as a function of time is characteristic of the attractor, and of the bifurcation routes 
through which attractors are created: intermittent dynamics leads to very distinctive distributions of local Lyapunov 
exponents [po[ . 

Our present results further underscore the utility of finite-time Lyapunov exponents in describing the local structure 
of dynamical attractors |2|]^,|20t in general. For the case of hyperbolic attractors, the theory connecting these to 
the invariant measure is well-developed. The present paper is part of a preliminary step towards understanding the 
connection between an underlying fractal structure and globally nonchaotic dynamics on strange nonchaotic attractors. 
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Fig. 1. Phase diagram for the quasiperiodically forced logistic map (schematic) [18], correspond- 
ing to the period 3 window of the unforced case. In order to obtain this diagram, A is calculated 
in a 100 x 100 grid. The shaded region (S) along the boundary (the A = contour) shows the 
region of SNAs. Intermittent SNAs are found on the edge of the Ci region marked I, while the right 
boundary, denoted C'l has fractalized SNAs. The left boundary between the periodic region and Ci 
does not show any SNA. The boundaries of the attractors can be interwoven in complicated manner 
(especially along C'l). 




Fig. 3. As in Fig. 2, but for N = 1000. The parameters for the four systems are the same in Fig 2. 
(a)-(d). A gaussian is fit at the maximum of the density, and is shown here as a dotted line. 




Fig. 4. Variation of F^(N) for the different SNAs, the parameters being identical to those in Figs. 
2. and 3., the fractalized route (□), the HH route (V), the blowout bifurcation route (x), and the 
intermittency route (o). 




Fig. 5. Variation of the (a) variance, u'^ and (b) the coefficient of skewness (71) with N for the 
different SNAs. The symbols and parameter values are the same as in Fig. 4. 



